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ABSTRACT 

We present evidence that the star-forming region NGC 346/N66 in the Small Magellanic Cloud is 
the product of hierarchical star formation, probably from more than one star formation event. We 
investigate the spatial distribution and clustering behavior of the pre-main sequence (PMS) stellar 
population in the region, using data obtained with Hubble Space Telescope's Advanced Camera for 
Surveys. By applying the nearest neighbor and minimum spanning tree methods on the rich sample of 
PMS stars previously discovered in the region we identify ten individual PMS clusters in the area and 
quantify their structures. The clusters show a wide range of morphologies from hierarchical multi- 
peak configurations to centrally condensed clusters. However, only about 40 per cent of the PMS 
stars belong to the identified clusters. The central association NGC 346 is identified as the largest 
stellar concentration, which cannot be resolved into subclusters. Several PMS clusters are aligned 
along filaments of higher stellar density pointing away from the central part of the region. The PMS 
density peaks in the association coincide with the peaks of [O III] and 8 /im emission. While more 
massive stars seem to be concentrated in the central association when considering the entire area, we 
find no evidence for mass segregation within the system itself. 

Subject headings: stars: formation — stars: pre-main sequence — galaxies: star clusters — Magellanic 
Clouds — Methods: statistical 



1. INTRODUCTION 

Star formation is generally thought to be clumped 
into a hierarchy of structures, from small multiple 
systems to giant stel l ar co mplexes and beyond (e.g. 
lEfremov fe Elmegreenl Il998h . and the interaction be- 
tween gravity and supersonic turbulent motions has been 
proposed as a major p layer in the formation of hierar- 
chical structures dMac Low fc Klessenl 120041: lElmegreen 
20061: iBallesteros-Paredes et al.ll2007t iMcKee fc Ostriker 
2007). The interstellar medium has also be en observed 
to be arranged in a hier archical structure (|Scalol Il985t 
IVazquez-Semad cni 2004), sometimes described as frac- 



tal, from the scales of the largest giant molecular clouds 
(GMCs) down to individual cores and embedded clusters, 
which are sometimes hierarchical themselves. 

The velocity structure in such turbulent regions is sug- 
gested to be scale-free, meaning that the same physical 
processes and the same velocity-separation r elations oc- 
cur o ver a wide range of absolute scales (e.g. Elmcgrcen 
2006). According to this hypothesis, large clumps orig- 
inate from the large velocities that compress the gas on 
large scales, and small clumps result from the compres- 
sion of the gas characterized by small velocities on small 



Electronic address 
Electronic address 
Electronic address 



sschmeja@ita.uni-heidelberg.de 
dgoulicr@mpia. de] 
rklcsseniaita. um-heidelberg. de 



scales. Considering that stars are formed by the col- 
lapse of such hierarchical structured clumps, they are 
themselves hierarchically struct ured. Indeed, stars ar e 
mostly formed in groups (e.g. IStahler fc Pallal [2005), 
whic h can be seen as the b ottom parts of this hierarchy 
(e.g. lElmegreen et al.ll2000D . 

Naturally, the structure of such systems, i.e. the spa- 
tial distribution of their members, may hold important 
information on the mechanism and the initial conditions 
of their formation. However, a thorough investigation 
of the structural behavior of star formation in environ- 
ments different than our Galaxy was not possible, until 
recently, due to observational limitations. In this paper 
we present such an investigation for the brightest star- 
forming region in the Small Mag ellanic Cloud (SMC). 

LHA 115-N66 or in short N66 (|Henizd fl956') is a very 
active star-forming region in the SMC with a diame- 
ter of more than 9 a t a distance of about 60kpc (e.g. 
lLanev fc Stobi3fi99l. Th is H II region is also known as 
DEM S 103 (jDavies et aL|[l976T) or NGC 346 referring to 
the bright stellar association located in its center. Being 
the largest and most luminous H II region in the SMC, 
NGC 346/N66 has been mapped in several wavelengths. 
Observations of its area have been taken in the X- - 
rays with XMM-Newton and Chandra (jNaze et al.ll2002L 
l2004h . in the UV with the Inte rnational Ultraviolet Ex- 
plorer (|de Boer fc Savagelll980l ) and the Far- Ultraviolet 
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Fig. 1. — Color-composite image of the three partially overlap- 
ping ACS/WFC observed fields of NGC 346, covering a total area 
of 5' X 5'. North is up, and east is to the left. Image Credits: 
NASA, ESA and A. Nota (STScI/ESA). 



Spectroscopic Explorer (|Danforth et alj 20031) . in the 
CO(2 -l) line with the SEST telescope (|Contursi et all 
2000), and in the infrared wi th the Infrared Space Obser- 
vatory (jContursi et al.ll2000f) and more recently with the 
Spitz er Space Telescope ( Bolatto et al.|[2007l : lSimon et al.l 
l2007f ). The region has also been covered by observa- 
tions of the radio continuum (lYe et al.lll991t iReid et al 



2006) and of the 21-cm atomic line ( Stavelev-Smith et al 
1997f ) . Ma ps of the ionized gas have also been made 
in Ha (e.g. iKennicut^ 11988^ Coarer et aTlll993| ). and 
[O III] with NTT (e.g. lRubio et al.ll2000f ). 

N66 is characterized by weak molecular emission and 
is also deficient in H I. It hosts, though, a variety of 
discrete H2 emission peaks, which coincide with the 
main features of the ionized gas, and with compact 
embedded young clusters, where candidate young stel- 
lar objects (YSOs) have been identified. The central 
part of the whole H II region is characterized by an 
oblique bright emission region, extending from south- 
east to northwest, known as the bar of N66 (Fig. [1]). 
This area is dominated by a strong UV radiation field 
related to the central association NGC 346. Indeed this 
system hosts the largest sample of s pectroscopically con - 
firmed OB stars in the whole SMC (jMassev et"al"1ll989f h 
which have b een the subject of s e veral previous inves- 
tigations (e.g. iNiemela et al.1 119861: IWalborn et "all 120001: 



lEvans et al.ll2006t iHunter et aLll2008ft . NGC 346 ionizes 
the remaining molecular cloud to the southwest of the 
bar, triggering the development of a photo-dissociated 
region seen as a curved front eating into the cloud. 

Recent observations with the Advanced Camera for 
Surveys (ACS) onboard the Hubble Space Telescope 
(HST) provided substantial information on the stel- 
lar content of this magnificent star-forming region by 
revealing a plethora of low-mass pre-main sequence 
(PMS) stars. This discovery clearly indicates that 
NGC 346/N66 also hosts low-mass star formation. Our 



photometric data set delivers more than 98 000 stars 
in total, with almost 8 000 being low-mass PMS stars 
(|Gouliermis et all [2006). The subsequent investigation 
of their spatial distribution showed that these stars are 
not homogeneously distributed across the region, being 
grouped in discrete concentrations, apart from the as- 
sociation NGC 346. Most of these "subclusters" are 
found to coincide with bright Ha emission, and the 
H2 emission peaks (jHennekemper et al.1 12008| ). Low- 
mass stars with masses 2-3 M^ remain in their 
PMS phase for ~ 2 0-30 Myr (e.g. iStahler fc Pailal [20051 : 
iBriceno et al.l l2"b07). As a consequence the rich sample of 
low- and intermediate-mass PMS stars, found with ACS 
in NGC 346/N66, offers a unique opportunity to investi- 
gate clustered star formation within a time-scale of the 
order of ~ 10 Myr and length-scales of 1 to 100 pc in 
a detail never achieved before for a star-forming region 
outside our Galaxy. 

In this paper we analyze the spatial distribution of 
the young PMS population unveiled by HST/ ACS in the 
bright SMC star-forming region NGC 346 /N66 and com- 
pare the results to current theoretical scenarios of star 
formation. For this purpose, we make use of our deep 
ACS photometry of this region, and apply two comple- 
mentary statistical methods to identify all the compact 
PMS clusters and analyze their structures. The data 
set used here is described in Section [2[ the statistical 
methods are explained in Section [3] and the results of 
our structural analysis of the PMS stars is presented in 
Section 0J A discussion on the star formation process 
revealed from our analysis is given in Section [5] 

2. DESCRIPTION OF THE DATA 

The observations used in this study, taken with the 
Wide-Field Channel (WFC) of ACS and their pho- 
tomet ry are thoroughly described by iGouliermis et al.l 
(2006).' The data were collected within the HST GO 
Program 10248 and were retrieved from the HST Data 
Archive 1 . Three overlapping images were taken with 
ACS/WFC centered on the association NGC 346 in the 
filters F555W (= V), F8UW (= I) and F6587V (= Ha), 
covering an area of about 5' x 5' (~ 80pc x 80pc at the dis- 
tance of the SMC). A color-composite image of this area 
is show n in Fig. [H It covers the intermediate-age c luster 
BS 90 dBica fc Schmittj [19951 : iRochau et al.l [2001 . the 
association NGC 346 itself, and all know n components 
of th e bright nebula N66 (jHenizd 119561 : iDavies et al.l 
1976). Photometry was performed using the ACS mod- 
ule of the package DOLPHOT 2 , and after eliminating 
bad detections based on the quality parameters of the 
detected sources returned from DOLPHOT, more than 
98 000 stars down to V ~ 27.5 mag were included in the 
photometric catalog. This catalog is complete ( 50%) 
down to V ~ 27 mag. In the present study, which deals 
with both low- and high-mass young populations in the 
area, we make use of the upper main sequence (UMS) 
stars with 12 mag ^ V ^ 22 mag, and the red PMS 
stars with V 22 mag, as they have been selected by 

1 The HST Data Archive is accessible from ESQ at 
http : //archive . eso . org/cms/hubble-space-telescope-data 
and MAST at http://archive.stsci.edu/hst/ 

2 The ACS mode of DOLPHOT is a n adaptation of the photom- 
etry package HSTphot (Dolphin 2000). The latest version can be 
downloaded from http://purcell.as.arizona.edu/dolphot/ 
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iGouliermis et all (|2006[ ). 

3. STATISTICAL METHODS 

3.1. Nearest Neighbor Density 

Star clusters are usually identified as regions of a cer- 
tain overdensity with respect to the background stellar 
de nsity. The nearest n eighb or (NN) method, introduced 
by iCasertano fe Hut! (|1985| ) based on earlier work by 
von Hoernerl (119631. estimates the local source density 
Pj by measuring the distance from each object to its jth 
nearest neighbor: 



Pj 



(1) 



where Tj is the distance of a star to its jth nearest neigh- 
bor, S(rj) the surface area with the radius rj and m the 
average mass of the sources (to = 1 when considering 
number densities). The value of j is chosen depending 
on the sample size. The chosen j value is correlated with 
the sensitivity to the density fluctuations being mapped. 
A small j value increases the locality of the density mea- 
surements at the same time as increasing sensitivity to 
random density fluctuations, whereas a large j value will 
reduce that sensitivity at the cost of losing some locality. 
Through the use of Monte Carlo simulations B. Ferreira 
& E. A. Lada (2008, in preparation) find that a value of 
j = 20 is adequate to detect clusters with about 10 to 
1500 members. 

The positions of the cluster centers are de- 
fined as the density- weighted enhancement centers 
(|Casertano fc Hutlll985l ) 



(2) 



where Xi is the position vector of the ith cluster member 
and p^ the jth NN density around this object. 

Similarly, the density radius r^ is defined as the 
density- weighted average of the distance of each star from 
the density center: 



rj. 



p) 



(3) 



It corresponds to th e observational core radius 
(|Casertano fc Hutlll985D . 

We identify clusters as regions that have 20th NN den- 
sity values 3ct above the average background level. 

3.2. Minimum Spanning Tree and Q 

The second method makes use of a minimum spanning 
tree (MST), a construct from graph theory, which is de- 
fined as the unique set of straight lines ( "edges" ) connect- 
ing a given set of points without closed l oops, such that 
the sum of the edge lengths is a minimum (|Boruvkall926t 
iKruskall Il95l iPriml \l951\ \. From the MST we derive 
the m ean edge length £mst- ICartwrigh^fc Whitworthl 
(2004) introduced the parameter Q = £mst/s, which 
combines the normalized correlation length s, i.e. the 
mean distance between all stars, and the normalized 
mean edge length ^mst- The Q parameter permits to 
quantify the structure of a cluster and to distinguish be- 
tween clusters with a central density concentration and 
hierarchical clusters with possible fractal substructure. 



Large Q values (Q > 0.8) describe centrally condensed 
clusters having a volume density n(r) oc r~", while small 
Q values (Q < 0.8) indicate clusters with fractal sub- 
structure. Q is correlated with the radial density expo- 
nent a for Q > 0.8 and anticorrelated with the fractal 
dimension D for Q < 0.8. The dimensionless measure Q 
is independent of the number of objects and of the cluster 
area. A detailed description of the method, in particu- 
lar its implement ation and normalization u sed for this 
study, is given bv lSchmeia fc Klessenl (|2006l ). 

4. RESULTS 
4.1. The Large- Scale Structure 

Figures andEb show the 20th and 35th NN density 
maps of the PMS stars in the NGC 346 field. While 
the first is used to identify clusters, the second traces 
the larger-scale structure, in particular, apart from the 
bar of N66, it emphasizes a filament (arm- like feature) 
of enhanced stellar density exte nding to the north-east 
from the central association ( see IGouliermis et al.ll2008l ) 
and another one running in east-west direction at the 
southern edge of the observed area. 

Figured: shows the 20th NN density of the UMS stars 
brighter than my = 18 in the same field. While the PMS 
stars show a highly inhomogeneous distribution with sev- 
eral readily identifiable clusters, the UMS stars show a 
distribution with only one central peak, which roughly 
coincides with the maximum NN density of the PMS 
stars. 

Figure [3] shows an image of NGC 346 in the [O III] line 
(501.1 nm) obtained with the ESO NTT overlaid with 
the PMS NN density contours. In the central associa- 
tion a high density of both gas and PMS stars is found, 
but the smaller clusters north and south of the center do 
not coincide with a higher gas density. The small clusters 
in the south appear to be aligned offset, but parallel to 
the gaseous filament. When comparing the PMS stellar 
distribution to the 8 /im observation of the region (Fig- 
ure |4j a remarkable correlation of the NN density peaks 
with the 8 ^m emission peaks is found for the central 
association and the clusters north of it. 

4.2. Identified Clusters 

Using the 20th NN density and the 3cr limit, ten PMS 
clusters are detected in NGC 346, labelled A to J from 
north to south (Fig. [2^,). 

The central association, Cluster F, is by far the largest 
cluster with 2349 PMS stars. Clusters A, B, D and E 
are parts of the filament pointing from the main associ- 
ation towards the north-east more clearly visible in the 
35th NN density map (Fig. Et>). Cluster C is a rather 
evolved open clu ster with a PMS populatio n at the age 
of ~ 5 — 15Myr (|Hennekemper et al.ll2008T ). and it may 
be not related to the much younger other clusters. While 
Clusters B and D show two clearly visible separate den- 
sity peaks, they are picked up as single clusters using 
the 3cr limit in the NN density maps. However, when 
using a 4cr threshold, cluster D breaks up into two clus- 
ters (Dl and D2), while of cluster B only the western 
peak re mains (B2). Clusters C, D, and E are also iden- 
tified bv iHennekemper et al.l (|2008l ) as their Clusters 1, 
2 and 3. Their Clusters 4 and 5 are part of the main 
association in the NN method. Clusters C, D, and E 
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Fig. 2. — NN density maps of NGC 346 (shown in a logarithmic 
scale): (a) 20th NN density of the PMS stars. The black lines 
indicate the cluster boundaries, defined as 3<r above the average 
background density. Also shown are the positions of the known O 
(diamonds) and B stars (triangles), (b) 35th NN density of the 
PMS stars, emphasizing the large-scale structure, (c) 20th NN 
density of the UMS stars brighter than my = 18. 
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Fig. 3.— NTT [O III] (501.1 nm) observation of NGC 346 over- 
laid with the PMS 20th NN density contours; the cluster boundary 
contour is shown by a thick line. 




Fig. 4.— Spitzer IRAC Channel 4 (8 (im) observation of NGC 346 
overlaid with the PMS 20th NN density contours (analogous to 
Fig.©. 

correspond to the s ubclusters Sc 16, Sc 14 + Sc 15, and 
Sc 13 identified by iSabbi et all (|2007t ). while the cen- 
tral association is broken into eleven subclusters (Sc 1 to 
Sc 11) by these authors. It is not possible to reproduce 
the number and l ocations of the central subclusters of 
ISabbi et all (|2007l ) using the NN method even with dif- 
ferent j values (5 < j < 50) and different threshold levels 
(2a < Pj < 6a). 

Table [T] lists the identified PMS clusters with the posi- 
tions of their 20th NN density centers (columns 2 and 3), 
the number of cluster members (column 4), their density 
radius (Eq. [3]) in pc (assuming a distance of 62kpc; col- 
umn 5), maximum 20th NN density (column 6), and Q 
(column 7). 
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TABLE 1 

Detected PMS clusters 



Cluster 
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( r . n -2\ 
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B 
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8.38 


0.59 


C 


14.87806 


-72.15229 


90 


1.79 


6.11 


0.74 


D 


14.81802 


-72.15552 


257 


2.92 


36.35 


0.65 


E 


14.79577 


-72.16745 


47 


1.21 


5.34 


0.92 


F 


14.76906 


-72.17724 


2349 


5.26 


40.02 


0.68 


G 


14.80844 


-72.20081 


139 


3.14 


3.45 


0.70 


H 


14.63330 


-72.20700 


18 


0.70 


2.49 


0.86 


I 


14.68797 


-72.21261 


71 


1.10 


13.22 


0.80 


J 


14.71087 


-72.21309 


26 


1.31 


2.83 


0.78 



Fig. 5. — Q versus density radius (upper panel) and number of 
stars (lower panel) for the identified PMS clusters (cluster F not 
shown). The dashed line at Q = 0.8 indicates the division between 
hierarchical and centrally condensed clusters. 

4.3. Clustering Parameters 

The detected clusters show a wide variety in the Q pa- 
rameter, ranging from strongly centrally condensed clus- 
ters with Q > 0.9 (clusters A, E) to highly hierarchical 
clusters with multiple density peaks and values Q < 0.7 
(clusters B, D). Also the central association (F) with its 
several density peaks shows, as expected, a hierarchical 
distribution (Q = 0.68). However, this value has to be 
taken with care, since due to a degeneracy of the Q pa- 
rameter, highly elongated stru ctures like this one sy stem- 
atically show lower Q values (jBastian et al.ll2008l ) . Ap - 
plying the correction suggested by iBastian et alJ (|2008l ) 
changes the value to Q s» 0.74, still in the hierarchical 
regime. To avoid these uncertainties, and to compare the 
values for PMS and UMS stars, we use a quadrilateral 
area around the central region with 14? 73 < a < 14? 83 
and -72?19 < 6 < -72?16 (shown as a white box in 



Fig. [2J:). Both PMS and UMS stars in this region are 
equally centrally concentrated (Q — 0.81). As far as the 
entire NGC 346 area is concerned, the UMS stars show 
the same Q value, while the PMS stars show a hierarchi- 
cal distribution (Q — 0.73). 

The subclusters of multi-peak clusters can be centrally 
concentrated or again hierarchical. The two subclusters 
Dl and D2 show Q values of 0.80 and 0.90, respectively, 
while B2 has a value of Q = 0.71. 

The Q parameters show a moderate correlation with 
the size of the cluster measured by and the number of 
stars, illustrated in Fig. This is in remarkable contra- 
diction to the embedded clusters found in nearby Galac- 
tic molecular clouds, where smaller cluste rs tend to have 
Q va lues well in the hierarchical regime (|Schmeia et al.l 
2008). This may be due to the fact that we are only 
tracing the more evolved PMS stars in NGC 346 that 
probably show a more centrally condensed structure than 
Class 0/1 proto stars, which are expe cted to be hierarchi- 
cally clustered (jSchmeja et al.l 12008? ) . 

4.4. Q versus Magnitude 

We investigate the variation of Q as a function of mag- 
nitude (and thus, mass) for the PMS and UMS sample. 
This is done for the entire area as well as for the cen- 
tral region. Q was calculated for sources brighter than 
a certain magnitude in steps of 0.5 mag as long as the 
sample contained at least 40 objects. The resulting Q 
versus magnitude plots for PMS stars and UMS stars 
are shown in Fig. [6] and respectively. Both PMS and 
UMS stars show a decline of Q with fainter magnitudes 
when looking at the entire area, indicating that the mas- 
sive stars are more centrally condensed than lower-mass 
objects. However, both samples show no correlation of 
Q with magnitude when considering the central region. 

4.5. OB Stars 

In Fig. [2k the positions of t he spectroscop i cally veri- 
fied OB stars in NGC 346 by iMassev et ail (|1989f) are 
shown overlaid on the 20th NN density map. While it is 
apparent that the OB stars are centrally concentrated in 
the association, it is interesting to note that they are not 
located in the densest (yellow-red) regions of PMS stars. 
This indicates that the density peaks of PMS stars ob- 
served in Cluster F may be an observational bias, due to 
the fact that the OB stars outshine the faint PMS stars in 
their vicinity. This provides further support to our iden- 
tification of the association as a single, not sub-clustered 
system from both observed UMS and PMS populations. 
As far as the clustering behavior of the OB stars alone 
is concerned, they do not seem to form any central con- 
densation. Indeed, when considering the 27 OB stars 
lying in the association, they show a distribution in the 
hierarchical range with Q = 0.77. 

5. DISCUSSION AND CONCLUSIONS 

5.1. Theoretical Concepts 

With the advancement of infrared detectors and the 
ability to do wide-area surveys, it has been recognized 
that a large fraction of all stars in the Milky Way form 
in clusters and aggregates of various size and mass scales 
and that isolated or widely distribu ted star formation i s 
the exception rather than the rule (|Lada fc L ada 2003). 
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Fig. 6. — Q versus limiting magnitude for the PMS stars in the 
entire area (upper panel) and in the central region (lower panel). 



Fig. 7. — Q versus limiting magnitude for the UMS stars in the 
entire area (upper panel) and in the central region (lower panel). 



It is therefore highly interesting to study the clustering 
behavior in neighboring galaxies as well. The complex 
hierarchical structure of molecular clouds provides a nat- 
ural explanation for clustered star formation. Molecular 
clouds vary enormously in size and mass. In small, low- 
density, clouds stars form with low efficiency, more or 
less in isolation or scattered around in small groups of 
up to a few dozen members. Denser and more massive 
clouds may build up stars in associations and clusters of 
a few hundred members. This appears to be the most 
common mode of star formation, at least in the solar 
neighborhood (jAdams fe Mv crs 2001). The formation of 
dense rich clusters with thousands of stars or more is 
rare. The cluster NGC 346, therefore, is a very extreme 
case of clustered star formation. 

Star formation is governed by the complex interplay 
between gravitational attraction in molecular clouds and 
any opposing processes such as supersonic turbulence, 
thermal pressure and magnetic fields. 

There are two main competing models that describe 
the evolution of the cloud cores. It was proposed in 
the 1980's that cores in low-mass star-forming regions 
evolve quasi-sta tically in magnetically subcritical clouds 
(|Shu et al.lll987t ). Gravita tional contraction is mediated 
by ambipolar diff usion (jMouschovias fc Spitzeri 119761 : 
lMouschoviaslll991h causing a redistribution of magnetic 
flux until the inner regions of the core become supercriti- 
cal and go into dynamical collapse. Measurements of the 
Zeeman splitting of molecular lines in nearby cloud cores, 
however, indicate magnetic field strengths that fall below 
the critical value, in some cases only by a small margin, 



in m a ny cases however, by factors of many (ICrutcherl 
fl99l iBourke et all l200lt iCrutcher et all l2008h . Mat^ 
netic fields therefore seem not to be the dominant agent 
regulating stellar birth. 

This has lead to the suggestion that the supersonic 
turbulence ubiquitously observed in molecular clouds 
is the main physical proces s governing star formation 
dMac Low fc Klessenl l2004t jBallesteros-Paredes et al.1 
I2007t iMcKee fc Ostrikerll2007t ). Supersonic turbulence 
plays a dual role. On large scales it can support clouds 
against contraction, however, on small scales it can 
provoke localized collapse. Turbulence establishes a 
complex network of interacting shocks, where dense cores 
form at the stagnation points of convergent flows. The 
density can be large enough for gravitational collapse to 
set in. However, the fluctuations in turbulent velocity 
fields are highly transient. The random flow that creates 
local density enhancements can disperse them again. 
For local collapse to actually result in the formation of 
stars, high density fluctuations must collapse on time 
scales shorter than the typical time interval between two 
successive shock passages. Only then are they able to 
'decouple' from the ambient flow and survive subsequent 
shock interactions. The shorter the time between shock 
passages, the less likely these fluctuations are to survive. 
Hence, the timescale and efficiency of protostellar core 
formation depend strongly on t he wavelengt h and 
strength of the driving source dKlessen etaD 120001 : 
iHeitsch et aIll200lHVazquez-Semadeni et al.ll2003| ). and 
accretion hist ories of individual protostars are stron gly 
time varying (jKlesseiJl200l iSchmeia fc Klessenll200l . 
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Interstellar turbulence is observed to be dominated by 
large-scale veloc it y mod e s (iMac Low fc O sscnkopf il200d ; 
lOssenkopf et all feoOlt lOssenkopf fc Mac Low! l2002f ). 
This implies it is very efficient in sweeping up molec- 
ular cloud material, creating massive coherent struc- 
tures which can go into large-scale gravitational col- 
lapse. In the extreme case, dense clusters with many 
thousand stars build up on one or two cross ing times 
(jHartmann et al.ll200H fElmegreenll2000l I2007D . Numer- 
ical simulations of gravoturbulent cloud fragmentation 
indicate that in such cases star clusters would build up 
in a hierarchical fashion with sub-clusters which evolve 
dynamically and merge toget her to build up the final , 
centrally-condensed cluster fKjessen fc Burkertl 120001: 
Clarke et al"1l200d: IBonnell et all 120031 : IClark fc Bonnelll 
20051 IBonnell fc Bate! I2006D , with stellar feedbac k 
strongly influencing this p rocess (|Krumholz et al.|[2007T ). 
iSchmeia fc Klesser]()2006l ) investigated the temporal evo- 
lution of a large set of numerical simulations of gravotur- 
bulent fragmentation and showed that the Q parameter 
of all resulting clusters increases with time. This is con- 
firmed by observations of embedded clusters in nearby 
molecular clouds, where the younger Class 0/1 sources 
show substantially lower Q values than the more evolved 
Class 2/3 objects (jSchmeia et al.ll2008ft . 

5.2. The Structure of NGC 346/N66 

Within our investigation of the spatial distribution of 
PMS stars in the NGC 346/N66 region using the near- 
est neighbor and minimum spanning tree methods an in- 
creased NN density of PMS stars is found in the central 
association and in two filamentary structures, one ex- 
tending to the north-east from the center of N66, and an- 
other at the southern edge of the observed area, extend- 
ing from east to west. We identify ten individual clusters 
of PMS stars in the region, with the central association 
being the largest. While several PMS density peaks are 
found within the association NGC 346, the system could 
not be resolved into individual subclusters. Our methods 
show evidence of centrally concentrated UMS stars and 
PMS stars in a hierarchical distribution. The PMS den- 
sity peaks of the association coincide with peaks of both 
[O III] and 8 /im emission, indicating that the system is 
in an early evolutionary stage. 

The identified PMS clusters in the whole area of 
NGC 346/N66 show a wide range of morphologies from 
hierarchical multi-peak configurations to centrally con- 
densed clusters. Only about 40% of the PMS stars in 
the investigated area belong to the ten identified clusters. 
The remaining 60% of PMS stars is found concentrated 
in small (statistically insignificant) groups and mostly 
dispersed over the entire area. This fraction of clustered 
PMS stars is comparable to that found in Galactic star- 
forming regions. In their SCUBA surv ey of the Perseus 
molecular cloud lHatchell et alJ (|2005l ) found that 40 - 
60% of the protostars and prestellar cores are located 
in small clusters (< 50 Mq) and isolated objects. In- 
vestigations in other Galactic star-forming regions like 
Ophiuchus, Orion and Monoceros yield a percentage of 
distributed populati ons between 20 and 60% with a "typ - 
ical" value of 25% (|Carpenteri 1200(1 lAllen et al.ll2007n . 
However, there are uncertainties due to the incomplete- 
ness of the samples and the definition of the clusters. 
Moreover, the clustering properties themselves depend 



on the age of the population considered and on the 
masses of the stars. As a consequence, it is very diffi- 
cult to compare and interpret these numbers. 

Based on the analysis of Q versus magnitude, both 
PMS and UMS stars show roughly the same degree of 
central concentration over the entire magnitude range. 
When we consider the entire area, Q is found to decrease 
with magnitude, providing evidence that both PMS and 
UMS stars are centrally concentrated in respect to the 
whole observed region of NGC 346. This, especially for 
the bright stars, is expected considering that the central 
association contains the largest fraction of stars. How- 
ever, when we consider only the central region of the area, 
Q does not seem to depend strongly on magnitude, and 
therefore no evidence of mass segregation in the central 
association is found. 

Whether the clu sters in NGC 346 have formed coevally 
dSabbi et al.ll200l or in a sequential process ov er sev- 
eral Mvr (jMassev et all 119891 : iRubio et al.ll2000l ) is still 



under debate. As young clusters seem to evolve from 
a hierarchical configuration to a centrally concentrated 
one with time, the fact that the identified clusters in 
NGC 346 show a wide range in their structures with 
Q values ranging from 0.65 to 0.93, may indeed indi- 
cate that they are at different evolutionary stages. How- 
ever, the absolute Q values are not correlated with the 
cluster age, but probably depend on a variety of addi- 
tional factors like initial condition s or turbulent energy 
in the cloud (jSchmeja et al.l [2008h . and therefore this is 
no direct proof that the clusters have formed at different 
times. However, the rather high Q values of the southern 
clusters (H, I, J) around or above 0.8 agree with other age 
indicators. While the peaks of the central concentration 
and the clusters to the north of it coincide with the peaks 
of dust emission (indicated by the 8 /im observations) and 
gas (traced by [O III]), the southern clusters appear to be 
associated with less amounts of gas and dust, suggesting 
an older age of these clusters. iHennekemper et all (|2008l ) 
find that the northern clusters comprise a younger PMS 
population (not older than 2.5 My r ) than the other clus- 
ters. Recently, iGouliermis et all (|2008h presented evi- 
dence that these clusters are the product of triggered star 
formation. In combination, all these indications strongly 
suggest that the identified PMS clusters are not the prod- 
uct of a single star formation event. 

In conclusion our investigation of the spatial distribu- 
tion of PMS stars in the region of NGC 346/N66 using 
the nearest neighbor and minimum spanning tree meth- 
ods showed that this region appears to be the prod- 
uct of hierarchical star formation, where the individ- 
ual clusters show different structures according to their 
size, evolutionary stage, position in the cloud and pos- 
sible additional internal or environmental factors. The 
overall structure and statistical properties of the re- 
gion are consistent with numerical c alculations of gravo- 
turbulent cloud fragmentation (e.g. iKlessen et al.1 120001 : 
IClark et al J 12005 IBonnell fc BateH2006h . NGC 346 /N66 
is a unique case, where the complexity of star formation 
can be observed in its entire length-scale. 
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